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ABSTRACT 

In light of recent results from the LHC, MEG and neutrino experiments, we revisit the issue 
of charged lepton flavour violation (LFV), which arises naturally in supersymmetric theories 
with massive neutrinos via flavour-violating soft supersymmetry-breaking masses for sleptons 
induced by radiative corrections. We link our results to the expectations for neutrino mixing 
angles, with particular focus on #13 in various Abelian flavour texture models, showing how 
LFV observables may distinguish between different possibilities. We focus on the radiative 
decays \i — > and r — > /i7, and on detection prospects at the LHC and a linear collider (LC). 
We use supersymmetric parameters favoured by — 2 and cosmological considerations that 
are consistent with LHC searches for supersymmetry and the Higgs mass. We find a class of 
scenarios where the LHC may be sensitive to LFV sparticle decays and LFV processes could be 
detectable at a LC with centre-of-mass energy above 1 TeV, whereas LFV lepton decays may 
be suppressed by cancellations in the decay amplitudes. 
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1 Introduction 



In recent years, the existence of neutrino masses and oscillations with near-maximal — u T and 
large u e — > mixing has been established by extensive input from atmospheric [I], solar [2] 
and long-baseline reactor [3] and accelerator [H [5] neutrino experiments. Initial input on the 
possible range of the #13 [6], [7] was provided by the T2K [8] and MINOS collaborations [9], and 
definitive evidence for a non-zero value of #13 has been provided by the reactor experiments 
Daya Bay [ID] and RENO [TT], and very recently also by Double Chooz [T2] . 

A natural expectation in theories with massive neutrinos is charged-lepton-flavour violation 
(LFV), which is enhanced in supersymmetric theories via the renormalization of soft supersymm- 
etry-breaking parameters. The link between neutrino oscillations and violations of the individ- 
ual lepton numbers L e>fi)T raises the prospect of observing processes such as fi — > ej, fi — > 3e, 
t — > fi'-f and fi — > e conversion on heavy nuclei |13j . 

The present experimental upper limits on the most interesting of these processes, summarised 
below, already constrain significantly the parameter spaces of theoretical models. 

BR(fx^e-f) < 2.4 x KT 12 [H], (1) 
BR(t->w) < 4.4 xlO" 8 [T5], (2) 
BR(t erj) < 3.3 x 10" 8 [15]. (3) 

The strongest constraint on radiative decays is the MEG upper limit on BR(/i — > ej) [14J. 

Within the supersymmetric framework, one should also keep in mind other possibilities for 
observing LFV processes, such as slepton pair production at a Linear Collider (LC) [161 El 
HH [191 EQl EU E2] and signals at the LHC [23l El E3 M, W\ EH [29], particularly in X 2 
X + e ± /i T X2 — > X + Ai ± T =F decays (here x is the lightest neutralino, assumed to be the lightest 
supersymmetric particle (LSP), and xi is the second-lightest neutralino). These decays could 
provide search prospects that are complementary to direct searches for the flavour-violating 
decays of charged leptons, particularly for heavy superparticle spectra. 

In this paper we re-evaluate the prospects for observable charged LFV, based on updated 
knowledge of neutrino mass and mixing parameters, including the recent measurement of #13. 
We work within the framework of the most natural mechanism for obtaining light neutrino 
masses, namely the see-saw mechanism [30], in which an effective Majorana mass matrix for 
light neutrinos, m e ff = ■ (M^)~ l ■ m® T , arises from Dirac neutrino masses of the same 
order as the charged-lepton and quark masses, and heavy Majorana masses Mjv- 

In supersymmetric theories, the neutrino Dirac couplings Y v renormalise the soft supersymmetry- 
breaking sneutrino and slepton masses, generating LFV naturally [21]. Even if the soft scalar 
masses were universal at the unification scale, quantum corrections between the GUT scale and 
low energies would modify this structure via renormalization-group running, which generates 
off-diagonal contributions. This effect is particularly interesting in see-saw models, where in 
general the Dirac neutrino Yukawa couplings cannot be diagonalized simultaneously with the 
charged-lepton and slepton mass matrices [31 J. Given the large mixing of the corresponding neu- 
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trino species, charged-lepton-flavour violation may occur at enhanced rates in supersymmetric 
extensions of the standard model, giving rise to observable LFV signals [131 E2l E31 EH l35| [36]. 

We analyse this possibility within the constrained minimal supersymmetric standard model 
(CMSSM) with universal scalar, gaugino masses and trilinear terms at the GUT scale, mo, 
Mi/2 and A , respectively, using textures for mass matrices that are inspired by GUT models 
with Abelian flavour symmetries [371 138] . These textures reproduce naturally the observed 
fermion mass hierarchies and mixing angles. Despite their phenomenological appeal, however, 
there are ambiguities and limitations due to the fact that the entries in the mass matrices are 
determined only up to 0(1) numerical factors. In a recent paper analysing textures derived by 
an SU(5) Froggatt-Nielsen mechanism, where all coefficients were treated as random variables 
in a Monte-Carlo scan of the multi-dimensional parameter space, it was shown that values of #13 
matching the recent data are statistically favoured [39]. Additional work in the same context 
may also be found in Ref. [40J. 

In this work, we extend these studies, showing that the neutrino parameters favoured by the 
data offer good detection prospects for detection of charged LFV in several channels. We 
study the complementarity between different radiative decays, deriving specific information on 
the allowed correlations between the decay rates. In all cases, we focus on sparticle spectra 
and values of tan/3 that are compatible with searches at the LHC, and also compatible with 
cosmological relic density considerations [H]. As an output of our analysis, we discuss detection 
prospects at the LHC and a Linear Collider (LC), taking into account the constraints that the 
recent LHC supersymmetry searches and the measurement of the Higgs mass impose on the 
observability of LFV processes. As we discuss later in more detail, one of the scenarios favoured 
by these constraints has relatively large right-handed sparticle mixing, and hence may lead to 
observable LFV in cascade decays to sleptons at the LHC [26], as well as at a LC with a 
centre-of-mass energy above 1 TeV. 

2 Neutrino Mass Textures inspired by SU(5) 

Over the recent years, a plethora of textures have been proposed to explain the data on neutrino 
masses and mixing. The new data on 6*13 provide additional constraints, excluding certain 
possibilities and constraining others. In our work, rather than reviewing the vast literature on 
the subject, we choose a representative model that fits the fermion data and is well-motivated 
on theoretical grounds. Nevertheless, we try to keep the results as generic as possible, placing 
emphasis on the links between physical observables. We also keep in mind that several a priori 
different theoretical models may converge to similar phenomenology, since they are matched to 
the same data. 

The example we choose is provided by a SU(5) GUT combined with family symmetries [371138]. 
The mass matrices are constructed by looking at the field content of the SU(5) representations, 
namely: three families of (Q,u c ,e c )i G 10 in representations, three families of (L,d c )i e 5 rep- 
resentations, and heavy right-handed neutrinos in singlet representations. This model therefore 
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has the following properties: (i) the up-quark mass matrix is symmetric, and (ii) the charged- 
lepton mass matrix is the transpose of the down-quark mass matrix, which relates the mixing 
of the left-handed leptons to that of the right-handed down-type quarks. Since the Cabibbo- 
Kobayashi-Maskawa (CKM) mixing in the quark sector is due to a mismatch between the 
mixing of the left-handed up- and down-type quarks, it is independent from the mixing in the 
lepton sector. In this way, it is easy to reconcile large atmospheric neutrino mixing with the 
observed small Vckm mixing angles. 

Within this framework, and following [57] ESI 1121 US], for example, the Yukawa matrices have 
the form 

e 6 e 5 e 3 \ ( e 4 e 3 e \ ( d 1 ^ e l 1±na l e^ 1±n ^ 

Y u oc ( e 5 e 4 e 2 , Y t oc Yj oc e 3 e 2 1 , Y v oc el" 1 ' | , (4) 

z 3 e 2 \ J \ e 3 e 2 1 / \ el ni l el n2 l e^ 

where Y u ^i, v stand for the Yukawa couplings of quarks, charged leptons and neutrinos respec- 
tively, and rii denote the U(l) charges of the heavy Majorana neutrinos. The heavy Majorana 
mass matrix is then given by 

~2|m| £ [ni+n 2 | glm+nal 

M N oc ( e 2 ^ £ \^+n 3 \ | _ ( 5 ) 

£ \ni+n 3 \ e \n 2 +n 3 \ £ 2|n 3 | 

There is no unique choice for the right-handed neutrino charges ni,n 2 ,n 3 , and several choices 
may be compatible with the low-energy neutrino data. We know, however, that the neutrino 
masses and mixing angles are related to the vl^l contributions in the effective neutrino mass 
matrix 

m eff « m u 1 - r rn v , (6) 
M N 

which, if calculated from the charges in (Hj), is of the form: 

E 2 

m eff oc ( e 1 1 | . (7) 




This form of m e ff is quite natural in the simplest see-saw models with a single expansion 
parameter and generic structures for the heavy and light Majorana mass matrices, due to 
cancellations that eliminate the dependences on the right-handed charges. Its predictions have 
been extensively analysed from a phenomenological point of view [H], and give a reasonable 
match to the data, provided there are no cancellations with large mixing in the charged lepton 
sector. Among other predictions, #13 turns out to be naturally of the correct order of magnitude. 
It is interesting to note also that, to lowest order in e, Y e Y^ has the same structure as m e ff, 
namely 

/ E 2 E E \ 

Y(Y} <xm eff oc £ 1 1 . (8) 
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The flavour mixing matrices are determined by the following diagonalizations of the Dirac and 
Majorana mass matrices: 



V e T (Y e Y/)V; = diag(y*,j£,j£), (9) 

V D T {Y v Yj)V* = diagfo^y^yj,), (10) 

U N T M N U N = diag(Mi, M 2 ,M 3 ), (11) 

U u T m eff U u = diag(m ul ,m U2 ,m U3 ) . (12) 

In terms of the above matrices, the Maki-Nakagawa-Sakata (MNS) matrix is given by 

U MNS = U = V}U„, (13) 

and can be parametrized as: 

U — V ■ diag(e-^ l/2 , e"^ 2/2 , 1) , (14) 

where 

(Cl2 c 13 S12C13 Si3e * 

— C23S12 — S23 s 13 c 12e l<5 C23C12 — S23Si3Si2e l<5 s 23 c 13 | , (15) 
S23 s 12 — C23Si3Ci2e l<5 — S23C12 — C23Si3Si2e l<5 C23C13 

and Cjj and Sy stand for cos% and sin^-, respectively. 



2.1 Predictions for Neutrino Observables 

Within the above framework of Abelian flavour symmetries there are ambiguities in the choices 
of coefficients, limited to a certain extent by requiring consistency with the experimental data. 
Here we proceed as follows: we take the expansion parameter to be e = 0.2 and multiply the 
entries of Xe, Y u and M N given by Eqs. (j3J) and (JSJ) by coefficients in the range ±[0.6,2] that 
we indicate by Uij and iVy for Ye, Y v and Mat, respectively. Given our ignorance of the 
CP-violating phase 5, we consider the case 5 = in the following: our results do not depend 
strongly on this choice. 

We analyse C(10 9 ) combinations of these coefficients, as follows: 

(i) We select the charged lepton Yukawa coupling matrix Y e so that the correct charged-lepton 
mass hierarchies are reproduced, namely 

^ ~ 0.06 — ~ 2.5 x HT 3 . (16) 
m T m T 

(ii) The coefficients and Ny of Y u and are then matched to a light neutrino mass 
matrix m e ff of the form ([7]), with entries that deviate by a factor between 1/2 and 2 from 
those in Eq. fl7J). 
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Figure 1: Correlations between the neutrino mixing angles before and after imposing the con- 
straints on the model coefficients as discussed in the text. The solid lines indicate the experi- 
mental bounds, and the small black crosses represent models satisfying all the constraints. The 
two large magenta crosses correspond to benchmark models that are discussed in the text and 
used for numerical calculations. 



(iii) We assume a neutrino mass hierarchy with a mass splitting of the order of the neutrino 
masses and m u , A ~ 0.05 eV. 

(iv) We require the following ranges for the mixing angles [TJ: 

0.27 < sin 2 <0.35, 
0.34 < sin 2 23 < 0.67, 
0.018 < sin 2 0i3 < 0.033. (17) 

The range on i3 is consistent with the values reported by both JTUJ [11] at the 3 — o 
level 0. 

(v) We make a further selection by requiring that the hierarchy of eigenvalues of Y(Y^ (which, 
as discussed above, has a similar structure to m e fj and Y u Yj) preserves the order of the 
gauge eigenstates. This simply reduces the density of solutions in the plots, since the 
predictions for the neutrino parameters are similar also for the other cases. 



In Fig. [T] we present the predictions for the neutrino mixing angles corresponding to the above 
criteria. We can see that, within this class of models, most solutions that reproduce the correct 
range of 0i2 and 023 also predict a neutrino mixing angle 0i3 that is compatible with the data 

1 Different choices of expansion parameters and coefficients iead to similar phenomenology. In ref. [39], an 
expansion parameter i] = 0.06 and coefficients in the range [-7=, vTo] are used. 
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Fit Ye Yv Mn 

/ e 4 1.5e 3 2e \ / e |1±ni1 e |1±n21 1.75e |1±n3 l \ / e 2 ^ £ N+n 2 | _ 1>3e |ni-hi 8 | \ 
1 I 0.5e 3 -1.9e 2 -0.5e s 1 ™ 21 1.25el" 3 l £ \m+m\ £ 2\n 2 \ £ \n 2 +n 3 \ J 

\ 2e 3 -1.9e 2 1.5 / \ l.3e^ -2eM ) \ -1.3sl ni+,l3 l £ l" 2 + n3 l - £ 2 M y 

sin 2 6*13=0.027 s j n 2 g 12 =0.30 sin 2 6* 23 =0.64 

2g|l±na| \ / g2|ni| g|m+»2| _ £ \ni+nz\ \ 

—0.5s' 113 ' I I £ l n l+ ,l 2l £ 2 M e |n 2 +n 3 | I 

1 25e'" 3 ' / \ — £ l ni + ,l3 l £ l»2+«3| _ £ 2|n 3 | / 

sin 2 6> 13 =0.019 sin 2 0i 2 =O.28 sin 2 fl 2 3=0.40 

Table 1: Indicative textures for Yi, Y u and M N , corresponding to the two crosses in Fig\J\ with 
the indicated values of the neutrino mixing angles. The are Abelian charges to be determined. 
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from [TOl [TT]. Furthermore, one can see that, within the narrow range of #13 indicated by 
the measurements of [10J E], the model predictions are in the mid-lower range of #12. On 
the other hand, in the case of 823, solutions are found in two disconnected regions in the plane 
(sin 2 # 13 ,sin 2 #23); & t either end of the range that is in agreement with a global statistical analysis 
of the data performed in However, we did not find any correlation between the preferred 
values of 823 and 9 12 in Fig. (TJ since the points in either of the two favoured regions of the right 
panel are distributed homogeneously in the left panel. 

We would like to point out that the fits in Fig. [U are independent of the charges n^, which 
affect Y v and Mn but not their combination in m e //. On the other hand, the choices of the 
do affect the rates for charged LFV processes, whose magnitudes provide important inputs on 
the right-handed sector of the neutrino mass matrices. Charged LFV processes are therefore 
powerful probes of parameters that cannot be constrained by lepton hierarchies or by better 
measurements of the neutrino parameters. We return to this issue later in the paper. 

As benchmarks for studying LFV in subsequent sections, we have identified two representative 
sets of textures from the many points that fit the neutrino data. These are indicated with crosses 
in Fig. [TJ Since we are working with an Abelian flavour symmetry, we focus on solutions where 
lepton mixing arises from large off-diagonal elements in Ve (in the basis of the GUT multiplets). 
In this case, the observed mixing is directly generated by the flavour symmetry - since it is 
correlated with the flavour charges of the charged leptons - and does not arise purely from the 
see-saw conditions, which often require more fine tuning. 

These benchmarks are chosen as representatives of the two different regions identified in Fig [TJ 
Fit 1 lies in the right region of the right panel of Fig [TJ with larger #23 , and Fit 2 lies in the 
left region of this panel, with smaller #23- Their parameters are listed in Table UJ For the 
reason discussed above, the mixing angles are insensitive to the values of the charges rii, which 
determine the levels of the LFV processes, as we see in the following sections. 
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3 Charged-Lepton-Flavour Violation 



We work in the constrained minimal supersymmetric standard model (CMSSM) with universal 
scalar, gaugino masses and trilinear terms at the GUT scale, m , M1/2 and A , respectively. The 
physical values of these masses are obtained by integrating the renormalization-group equations 
(RGEs) of the MSSM supplemented with right-handed neutrinos. Since the Dirac neutrino and 
charged-lepton Yukawa couplings cannot, in general, be diagonalized simultaneously, and since 
both lepton Yukawa couplings appear in the RGEs, the lepton Yukawa matrices and the slepton 
mass matrices at low energies cannot be diagonalized simultaneously, either. 

In the basis where the charged-lepton masses me are diagonal, the soft slepton-mass matrix 
acquires a correction that contains off-diagonal contributions from renormalization at scales 
below Mqut, of the form [32J: 

H « I^ 6m o + 2 ^o)^ln(^) , (18) 

where is the intermediate scale where the effective neutrino-mass operator is formed. The 
physical charged slepton masses are obtained by numerical diagonalization of the following 6x6 
matrix: 

«?i=(t? 5* ) ' (19) 

" L RL RR 



where all the entries are 3x3 matrices in flavour space. Using the basis where Ye is diagonal, 
it is convenient to write the 3x3 entries of ([19]) in the forms: 



<l = {mpf + 5m} + m}- l -{2M 2 w -Ml) C os2^ (20) 
<r = (m^) 2 + m, 2 - (Mf -M^) cos 2/3, (21) 



m 



RL 



(Af ag + 5A e -fi tan (3)m e , (22) 



m LR = m R L - ( 23 ) 



Here tan [5 is the standard ratio of the two MSSM Higgs vevs, (m~ iag ) 2 , (m~ iag ) 2 and v4^ iag denote 



the diagonal contributions to the corresponding matrices, obtained by numerical integration of 
the renormalization-group equations, and Sm'j and SAe denote the corrections with off-diagonal 
terms that appear because Y v and Ye cannot be diagonalized simultaneously. 

The full mass matrix for left- and right-handed sneutrinos has a 12 x 12 structure, given in 
terms of 3 x 3 Dirac, Majorana and sneutrino mass matrices. The effective 3x3 mass-squared 
matrix for the left-handed sneutrinos has the same form as the m\ L part (|23|) of the 6x6 
charged-slepton matrix ([191) . with the difference that now the Dirac masses are absent. In 
Ref. [23] it was shown that is sufficient to use 

ml = ( m f ag ) 2 + 5m) + -M 2 Z cos 2(3 . (24) 

1 2 
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The matrix responsible for LFV is Vlfv — V^Ve, and the slepton mass matrices contain off- 
diagonal terms generated by: 

m LL = ViFv( m LL)diagV L FV , (25) 

while the A-terms become: 

A £ = V? FV (A e ) dias . (26) 

Here (m| L )di a g and (A^^g are the soft terms resulting from the RGE running of the universal 
soft terms at the GUT scale. 

The matrix element of the electromagnetic-current operator between two distinct lepton mass 
eigenstates £i and £j is given in general by 

T A = (£i\( P - q)\JxMp)) 

= Ui{p- q) [m 3 ia X pq p (A M P L + A M P R ) + 

{q 2 lx ~ qxl • q) (A L E P L + A R E P R )] U] { P ) , (27) 

where q is the photon momentum. The coefficients Am and A E denote contributions from 
neutralino/charged-slepton and chargino/sneutrino exchanges. The amplitude of the LFV pro- 
cess is then proportional to T\€ X , where e A is the photon-polarization vector. The branching 
ratios (BRs) of the decays tj — » £i + 7 are calculated using the formulae of Ref . [32] . 

3.1 Predictions for Radiative Decays and LFV at LHC 

We now study the conditions under which the favoured ranges of lepton masses and of neu- 
trino mixing, can lead in a natural way to observable signatures for charged LFV. As already 
discussed, while the neutrino parameters are independent of the charges n i; this is not the case 
for LFV. As a result, LFV can provide a way to probe the right-handed neutrino sector, for 
which only limited information is available. 

In a previous section, we identified two representative benchmark fits suitable for studying the 
expected rates for charged LFV. As we saw previously, the level of charged LFV is determined 
by the product Vlfv — Vd ' anc ^ ^ nus by ^ ne charges which enter in the Dirac neutrino 
mixing matrix V^. However, in general, different choices of rij lead to different solutions for a 
large Vi in the GUT eigenstate basis, and it is possible to group these solutions in the following 
categories: 

• Solutions with a large Ve lead to large LFV, for all charge combinations that generate a 
small Vp. 

• For solutions with large V^, the mixing may add to or subtract from that of the leptons, 
depending on coefficients and phases. In the first case, large LFV is expected, whereas in 
the second case LFV is suppressed. However, it is difficult to find suppression in all three 
relevant mixing angles simultaneously. 
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(0 

{iii = 1, «2 = 0, ri3 = 0} 





(iv) 

{iii = 0, ri 2 = 1, 113 



/ 0.8(35 -0.385 -0.451 \ / 0.805 -0.385 -0.452 \ / -0.820 0.343 0.458 \ / 0.806 -0.401 -0.436 \ 
V L fv 0.182 0.885 -0.429 -0.064 0.700 -0.711 0.494 0.829 0.263 -0.437 -0.899 0.016 

\ 0.565 0.263 0.782 / \ 0.590 0.601 0.539 / \ -0.289 0.441 -0.849 / \ -0.399 0.178 -0.901 / 



Table 2: Values of Vlfv with e = 0.2 for Fit 2 of Tabled with large Vg. Fom left to right, we 
see the following features: (i) Vd is small, leading to large Vlfv, (H) an d (Hi) Vd is large and 
so is Vlfv, (iv) Vd is large and cancellations with Vt occur in the 2-3 sector, suppressing LFV. 
In the last two examples, the heavy Majorana mass matrix has an inverted hierarchy. 

These observations indicate that, even before looking at the influence of the charges rij or 
performing specific computations, LFV processes can be quite large for natural choices of 
textures. 

A survey of the dependences of the entries in Vlfv on the different right-handed neutrino 
charges is given in Table |2J For simplicity, we focus on Fit 2, noting that similar results hold 
for Fit 1. We see explicitly that different values of the charges result in different levels of 
LFV, even though the m e ff eigenvalues are unchanged, thanks to the see-saw conditions. 

In addition to the flavour-mixing effects, the rates for flavour non-conserving decays are sen- 
sitive to other physical quantities. In the supersymmetric GUT context explored here, they 
depend on the masses of the sparticles that mediate the flavour non-conserving processes. In 
the CMSSM studied here, we parametrize these masses in terms of the universal GUT-scale 
parameters m and Mi/ 2 , and use the RGEs of the MSSM to calculate the low-energy sparticle 
masses, taking into account low-energy threshold effects. Other relevant MSSM free parameters 
are the trilinear coupling A, the sign of the Higgs mixing parameter /z, and the value of tan/3. 
Here we choose values of the CMSSM input parameters that are consistent with cosmological 
relic density constraints, as well as LHC results. 

The recent LHC measurement of the Higgs mass H7] imposes severe constraints in the 
CMSSM parameter space. Specifically, a Higgs masses rrih ~ 125 GeV implies, in general, a 
relatively heavy sparticle spectrum that is consistent with the cosmological constraint on the 
neutralino relic density only in limited regions. A global analysis of the CMSSM parameter 
space was performed in [JT] , and yielded two almost equally good fits to the available data, one 
with relatively low sparticle masses and tan f3 ~ 16, and the other with larger sparticle masses 
and tan/3 ~ 45, corresponding to the following CMSSM sets of parameters: 



Both these choice of input parameters are compatible with measurements, within theoretical 
uncertainties jH]. 

Point (a) predicts a lighter sparticle spectrum, which generally favours larger LFV rates at 



(a) tan /3 = 16, m = 300 GeV, M 1/2 = 910 GeV, A = 1320 GeV , 
(6) tan/3 = 45, m = 1070 GeV, M 1/2 = 1890 GeV, A = 1020 GeV . 



(28) 
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Figure 2: Predictions for the rare LFV decays £{ — > iff as a junction of the heaviest right- 
handed neutrino mass Mat, for the benchmark point displayed in (28\) (a) (thick line), (b) (thin 
line) and using some of the neutrino mixing fits shown in Tabled The solid lines correspond 
to cases (iv) and dashed lines (Hi). The horizontal solid lines indicate the current experimental 
upper bounds. 

accelerators, since there is more phase space at the LHC and smaller energies are required for 
sparticle production at a LC. On the other hand, point (b) is indicative of models with a heavier 
spectrum, presenting greater experimental challenges. Regarding cosmological considerations, 
point (a) belongs to the region where the WMAP-favoured range of Q x h 2 is achieved via x ~ ^ 
coannihilationa, whereas point (b) lies in the funnel region where the neutralino LSP annihilates 
rapidly via direct-channel H/A poles. 

2 We note that in this region the x ~ t mass difference is very small, offering other experimental challenges 
and opportunities [48] [49] . 
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In Fig. [5] we show numerical predictions for the LFV branching ratios arising from the textures 
introduced above (which match the neutrino data) for the CMSSM parameters specified in 
(128]) (which are compatible with the LHC constraints). We show the effect of varying 
from 6 x 10 14 GeV down to 10 12 GeV, for Fit 2 of Tableland for two choices of right-handed 
neutrino charges, rii (0,1,0) and (2,0,1) respectively). We can see that the experimental upper 
bounds on BR(yU — > ej) can be reached for suitable values of the Yukawa structure, even with 
the heavy sparticle spectrum suggested by the higher-mass benchmark (6). 

It is clear from the three panels of Fig. [2] that the combined uncertainties from the unknown 
right-handed charges rii and the unknown model coefficients inevitably introduce significant 
uncertainties. On the other hand, the rates for radiative charged LFV exhibit different correla- 
tions in the different benchmark scenarios, suggesting that future observations could possibly 
provide important input into the selection of model parameters. Further information could be 
obtained by LFV studies at a Linear Collider (LC), as discussed in the next Section. In general, 
we find that: 



• The branching ratios are larger for the lower-mass scenario with tan/3 = 16, due to the 
lighter spectrum. 

• Large values of Aq may induce some cancellations in the branching ratios for radiative 
decays, suppressing them strongly. In such regions, colliders provide the best chances for 
probing LFV. We discuss some examples at the end of the next Sections. 

• A lighter spectrum implies smaller threshold energies and larger cross sections for slepton 
production at a LC. 



In previous work [26], we concluded that the bounds of BR(r — > 717) disfavour mixings that 
could lead to observable signals at the LHC in the LL sector, whereas significantly larger mixing 
were still possible in the RR sector. This occurs for a broad region of points like scenario (6), 
where the size of A is limited by electroweak symmetry breaking and rrih ~ 125 GeV is realized 
with a heavy spectrum with large tan/3. Point (a), on the other hand, has a light spectrum with 
large trilinear terms that are nevertheless compatible with electroweak symmetry breaking. As 
a result, LR flavour mixing in the charged slepton masses may induce non-trivial LFV in slepton 
decays. 

In what follows, we revisit the X2 — > X + T ± Ai =F decay mode for the benchmark points (a) and 
(6), in order to determine whether the relevant ratios 

Rr, = T( X 2 -> X + t± + /i T )/r( X2 -> X + t± + r T ) (29) 

can be of the order of 10%, as required for a signal that can be distinguished from the back- 
ground [24j|26]. We include all the contributing on-shell sfermion exchange diagrams, as given 
in EDI: 



BR(x 2 -> XtV) = [ BR (X2 -> £/x)BR(£ -> t X ) + BR(x 2 -> £t)BR(£ -> //*)] . (30) 
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Figure 3: The ratio defined in Eq. [29) for the CMS SM points (a) (thick line) and (b) (thin 
line), Eq. H21 using the mixing options (iv) (solid lines) and (Hi) (dashed lines) of Tabled 

We present in the left panel of Fig. [3] the ratio of the flavor- violating decay width T(x2 — > 
X + T± + H T ) to the CMSSM flavor-conserving decay width T(x2 -* X + r± + r=F )- We see that 
for point (a) predictions are within the reach of the LHC, whereas for point (b) the predictions 
are below the expected experimental sensitivity for LFV observation in neutralino decays at 
the LHC. 



4 LFV Observation at a Linear Collider 

In previous work [21], we studied the different LFV channels at a LC, focusing on the production 
of /i ± r ± pairs. In supersymmetric models where LFV is produced by lepton-slepton vertices, 
observable signatures may occur either directly, in slepton-pair production, or indirectly, via 
slepton production in cascade decays [TB] . If the flavour mixing is introduced in the LL slepton 
sector, as is the case for the models under consideration here, pair production dominate in 
slepton-pair production and LFV decays such as: 



At the CMSSM benchmark point introduced above, which favoured by current experimental 
data and cosmological observations, the channel mediated by charged sleptons clearly dominates 
over the sneutrino-pair production process, and may lead to a cross section of the order of 1 fb, 
the reference value used in [21], for a LFV signal of /i ± r ± pairs that can be distinguished from 
the background, according to the considerations in [19] . Here, we extend our previous results 
by considering the full structure of the Yukawa couplings, thus comparing the LFV production 




(31) 
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of charged leptons of all generations. In addition, we update the supersymmetric benchmark 
points used previously, choosing a new CMSSM benchmark scenario compatible with the LHC 
measurement of m h . 
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Figure 4: Values of the cross sections <r(e + e — > £ i — > £^£^ + 2%°) (£ a ^ £ b as indicated 
for each panel) as functions of for Mn = 10 14 GeV. The line styles as the same as those in 
Fig. 03 

We present in Fig. @] the expected cross sections as functions of a/s for Fits 1 and 2 of Table 
[Hand the same sets of charges as those in Fig. 2. For defmiteness, the right-handed neutrino 
mass scale is set to = 10 14 GeV. We observe that, by varying this scale in a similar range 
as in Fig. [2} the cross sections also change by approximately an order of magnitude: according 
to the see-saw conditions, larger values of imply larger Dirac neutrino Yukawa couplings, 
and thus stronger RGE running. The results depend sensitively on the sparticle spectrum. If 
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Figure 5: Upper panels: branching ratios for the CMS SM points (a) (thick line) and (b) (thin 
line), Eq. ( ffgj) . using case (ii) the neutrino mixing fits shown in Table IB The lower panels 
show the corresponding prediction for the LFV cross sections. 

we restrict ourselves to the benchmark CMSSM parameter choice discussed above, we would 
require y/s beyond 3 TeV in order to have significant LFV signals in slepton mixing. However, 
once above threshold, the uncertainties in the choices of right-handed neutrino charges and 
unknown model coefficients have similar impacts on the LC LFV cross sections to those on the 
radiative decays discussed in the previous section. 

In the case of point (a) we find cancellations in the branching ratios for some of the neutrino 
mixing fits. This happens because of the large value of A Q , which leads to significant cancella- 
tions among the different amplitudes that contribute to the LFV decay amplitudes [32 | |33 | IM]. 
Using fits (ii) and (iv) from Tableland the benchmark point (a), we find this type of cancel- 
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Figure 6: The ratio defined in Eq. [29) for the CMS SM points (a) (thick line) and (b) (thin 
line), Eq. ( El| ), using the mixing options (ii) of Tabled 

lation in BR(r — > fi'j) and BR(r — > ej) for a range of values of the right-handed neutrino scale 
Mjy, as shown in the upper panels of Fig. The lower panels show the corresponding predic- 
tions for LFV stau decays at a LC when Mjv = 2 x 10 13 GeV. We see that the cancellations 
in the branching ratios do not occur for the LFV LC signals. This implies that it is possible 
to observe slepton flavor oscillations at the LC, in cases where LFV would be undetectable in 
rare charged lepton decays. 

Fits (z) and (ii) in Table [2] correspond to cases with larger mixing in the matrix Vlfv- m these 
two cases, the choice of charges induces an inverse hierarchy in M^. In both cases, the large 
yl-terms and large 2-3 mixing make cancellations possible. We display results for case (i) but 
not for case (ii), which is quite similar. 

5 Conclusions 

Abelian flavour symmetries provide interesting possibilities for understanding the hierarchy 
of fermion masses and mixing. Despite uncertainties in the choice of 0(1) coefficients, they 
offer useful insight into physical observables, and provide specific predictions for the signals to 
be expected in various detection channels, which serve as diagnostic tools for discriminating 
between different models. 

In our work, we have explored these issues using updated experimental input from neutrino 
data, particularly recent measurements of #13, MEG and the LHC. We have revisited the 
signatures of charged LFV in theoretically-motivated scenarios, studying the correlations arising 
in CMSSM models with parameter values that are favoured by the LHC and cosmological 
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considerations. Because of their sensitivity to flavour symmetries and model parameters that 
are not constrained by the neutrino data, particularly those linked to the right-handed neutrino 
sector, LFV searches may become a powerful tool for distinguishing between different theoretical 
scenarios. 

Within an SU(5) GUT framework supplemented by an Abelian flavour symmetry, we first 
performed a scan of different fits to the neutrino data, paying attention to the the naturalness 
of the fit and avoiding artificial cancellations arising from specific choices of coefficients. This 
leads to a pattern of correlations between the neutrino mixing angles that are compatible with 
the global analysis of neutrino data in [4"5] . 

We then looked at the expectations for LFV processes in the above models, identifying the 
ranges of parameters where observable signatures are possible. In general, fits with similar 
predictions for the neutrino parameters may lead to different LFV predictions. However, the 
recent input on 13 , combined with the new MEG bound on \i — > as well as LHC data, do 
constrain the allowed structures. Further precision in the determination of neutrino parameters 
could lead to restrictions on the choices of model coefficients, but would not constrain the right- 
handed neutrino charges. However, new input in this respect could be provided by the rates for 
LFV processes, since their magnitudes are directly linked to these charges, unlike the neutrino 
mass and mixing parameters. 

In the cases we studied, it was possible to establish correlations between the expected rates 
for radiative LFV decays, LFV processes at the LHC and r/x pair production at a future LC, 
for different possibilities for the structure of the heavy Majorana neutrino masses. As we have 
shown, extending previous studies by including the full 3x3 neutrino mixing leads to interesting 
correlations between the LFV production of charged leptons of different generations at a LC. 
However, within the CMSSM framework studied, the absence of a supersymmetry signal in the 
LHC data and the discovery of a neutral Higgs weighing ~ 125 GeV imply that observation 
of slepton flavour violation at a LC will be possible for only energies beyond 1 TeV, in models 
with large Dirac neutrino couplings and a heavy Majorana neutrino scale M^. On the other 
hand, scenarios less constrained than the CMSSM that could fit with a lighter sparticle 
spectrum might predict observable LFV signals at smaller energies. 
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